MACDONALD'S IDENTITIES AND THE LARGE N LIMIT 
OF YM 2 ON THE CYLINDER 



STEVE ZELDITCH 



Abstract. The purpose of this paper is to determine the large N asymptotics of the free 
energy Fjy(a,U\A) of YM 2 (two-dimensional Yang Mills theory) with gauge group Gn — 
SU (N) on a cylinder where a is a so-called principal element of type p. Mathematically, 

F N (U U U 2 \A) = logH GN (A/2N, U u U 2 ) 

is the central heat kernel of Gn- We find that 

N 

F N (a Nl U N \A)~ — E(d9,da) 

where 3 is an explicit quadratic functional in the limit distribution da of eigenvalues of 
Uni which depends only on the integral geometry of SU(2). The factor of N appears to 
contradict some predictions in the physics literature on the large N limit of YM 2 on the 
cylinder (due to Gross-Maty tsin, Kazakov-Wynter others). 



1. Introduction 

The purpose of this note is to determine the large N asymptotics of many cases of the 
free energy -Fjv(^Cn Uc 2 \A) of YM 2 (two-dimensional Yang Mills theory) with gauge group 
SU(N) on a cylinder. Interest in this large N limit problem was raised around ten years 
ago in a series of physics papers by Douglas, Kazakov, Wynter, Gross, Matytsin and others 
DK. \KW\ IGM1| IUM2 . They predicted, on the basis of physical and formal mathematical 
arguments, that the large N limit of F N {Uc 1) Uc 2 \A) should be well-defined and related to the 
action along a path of the complex Burgers equation with boundary densities determined by 
the limiting eigenvalue densities of Uq 1 , Uc 2 ■ In certain cases they predict a phase transition 
between a weak coupling and a strong coupling regime. 

Our results give the asymptotics of Fjv(ajv, Un\A) where is a principal element of type 
p of SU(N) in the sense of jKoj (see Section |2J), and where the second sequence {Un} can 
be any sequence of elements of SU(N) which possesses limiting eigenvalue density da. The 
limiting eigenvalue density of is uniform measure dd on S 1 . Our asymptotics reveal some 
surprising results: 

• Fn(Uc 1 ,Uc 2 \A) ~ 4 ^(d9,da) for a certain (explicit) functional S, unlike the pre- 
dictions that it tend to a limit functional as N — > oo. The extra factor of N signals 
an error in some of the heuristic physics arguments, and (as emphasized to the author 
by M. Douglas and C. Vafa) calls into question whether the large N of YM 2 on the 
cylinder is actually a string theory as discussed in |GMlj IGM21 IKWj and elsewhere. 
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At least, it indicates that extra hypotheses on the matrix pairs Uc x , Uc 2 are necessary 
for the conjectured picture to be correct. 
• The principal asymptotic term is visibly analytic in A, i.e. it never exhibits a phase 
transition to leading order, even though da could be any probability measure on 
S 1 . This aspect of the results is consistent with the physics predictions, in that the 
boundary condition dQ puts the system in its strong coupling regime and it should 
therefore not exhibit a phase transition between weak and strong coupling [D] . 

The calculation presented here appears to be the first rigorous calculation of the large TV 
limit of the partition function for YM 2 on the cylinder. Hence, it is unclear at this time 
whether the anomalous factor of iV only occurs in the special case where Uc x = cln is a 
principal element of type p or whether it holds more generally in this context. 

We should emphasize that the other rigorous mathematical results (of which we are aware) 
have largely confirmed the large N limit picture developed in |DKt IGMH IGM2| iKWj and 
elsewhere. The Douglas-Kazakov phase transition of the genus partition function at A = it 2 
has been proved by Boutet de Monvel - Shcherbina in |BSj . A large deviations analysis of 
spherical integrals by Guionnet-Zeitouni in |GZj justifies some of the predictions by Matytsin 
[M] and Gross-Matytsin |GMlj on the asymptotics of characters Xr{U)- However, the results 
of |GZj pertain to the analytic continuation of xr to positive elements of GL(N, C) and it is 
possible that Matytsin's predictions fail when U = ajv- We intend to explore this question 
in a future article. 

To state our results, we introduce some notation. The partition function of YM 2 on a 
cylinder, with gauge group equal to G, is given by Migdal's formula (see [Wl IW2j ) : 

Za{U 1} U a \A) = Y^XR^XRme^ ^ (1) 
Red 

Here, A > is the area of the cylinder, and the sum runs over the irreps (irreducible 
representations) of G, with xr the character of R and with C 2 (R) equal to the eigenvalue 
of the Casimir A in the irrep R. Thus, Zq is the value at time t = ^ of the central heat 
kernel of G: 

H a {t,Ui,U 2 ) = $>*(£/i)x*(£/ 2 *)e-' C2 ^ , (2) 
Red 

i.e. the kernel of the heat operator acting on the space of central functions on SU(N). It 
is obtained from the usual heat kernel by averaging both variables over conjugacy classes. 
Since Zq(Ui, U 2 \A) is conjugacy invariant, one may assume that Ui, U 2 are diagonal and we 
write Uj = -D(e l6, i, . . . ,e 1 ^), where £ is the rank of G (i.e. the dimension of its maximal 
torus). The main quantity of interest is the free energy, defined by 

F G (U Cl ,U C2 \A) = j^lnZ G (U Cl ,Uc 2 \A)- (3) 

We now consider the large N limit of the partition function and free energy. The large 
N limit refers to an increasing sequence Gn of groups, e.g. the classical groups G^ = 
U(N), SU(N), SO(N), Spin(N). For the sake of simplicity we restrict attention to SU(N) 
and we abbreviate Fsu(n) by (etc.) The limit we are interested in is a pointwise limit of 
the central heat kernel, which obviously requires some discussion since the space on which 
the heat kernel is defined changes with N. 
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The large N limit of Fn is defined as follows: take a pair of sequences {UNj},UNj £ SU (N) 
(j = 1,2) of elements whose eigenvalue distributions 

da Nj :=—J2^ <3 ) G MS 1 ) 



lN k=i 



tend to a limit measures o r Here, £n = N — 1 denotes the rank of the group, and M^S 1 ) 
denotes the probability measures on the unit circle, and convergence is in the weak sense of 
measures. We will denote this situation by U^j —* <J r That is, 

, e N 

U Nj - <7j <=► -^5(e< J ) - a E MiS 1 ), (j = 1,2). (4) 



Conjecture 1.1. |GM1| (pages 8-9) Assume that U\ —> ai,U 2 — > cr 2 . Then 
F N (U 1 ,U 2 \A)^F(cr 1 (0),<r 2 {O)\A) = S(* 1 (e),<T2(e)\A) 

-\ J sl J sl (TxiB)^) log | sm^\d6d<f> - \ J sl J sl a 2 {9)a 2 {(f>) log | sin -^\d6d(f>, 
where the functional S is a solution of the Hamilton Jacobi equation 

3A = 2f ^ws^e? -T"M ] - 

It is often assumed in the physics articles that the limit measures have densities. The 
function S can be identified as the action along the solution of the boundary value problem 
for Hopf-Burgers equation 

at ~ J dx u ' 

%f(0,x)=a 1 , Zf{A,x)=a 2 . 

Under certain assumptions on the limit densities, it is further conjectured that a third 
order phase transition between the weak and strong coupling regimes should occur. In the 
case of the disc (the cylinder with U\ = I), Kazakov-Wynter and Gross-Matytsin have 
predicted a phase transition point at 

A=-^-, 9(61):- fda{e ' ] 



26(71-)' v ' J 9 — 9' 

Here, da is the limit distribution for U 2 and the integral is presumably to be understood 
as the Hilbert transform on S l (although it is written as the Hilbert transform on M). The 
weak coupling regime is characterized by a gap in the support of the limit density (i.e. an 
interval of S 1 on which do\ vanishes). 

The mathematical evidence for these conjectures is largely based on approximating the 
discrete sum (JTJ) over R e G by a continuous (and only partially defined) integral over a 
space of "densities of Young tableaux" , to which the saddle point method is applied (cf. 
§3). These arguments are not rigorous and also make some implicit assumptions (see the 
discussion after Theorem II. 3|) . 

The main observation of the present article is that there exists a simple, rigorous alternative 
for calculating the large N limit when one of the arguments is a so-called principal element 
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of type p. This alternative method is based on the use MacDonald's identities (in Kostant's 
formulation) to factor the partition function as a product over positive roots. Background 
on MacDonald's identities will be provided in 32 (see also [MAC! iKol IF1 IPS]). 

The first result one obtains this way is a limit formula in which the time variable in the 
partition function is not rescaled. We use the notation da * da(x) for the measure defined 
by 

f{e ix )da*da{x):= [ [ f(e i{x - y) )da(x)da(y). (5) 
Js 1 Js 1 

Thus, * denotes convolution of measures and da is short for da(x). For notational simplicity 
we sometimes denote e tx G 5 1 more simply by x. 



Theorem 1.2. Let e su(iV) be diagonal matrices with entries 9f , and assume da 
± E$=i <Ke^) -> <J. Then, as N -> oo, 



N 



j^log H S u(N)(t,a N ,e kN ) -> -\ log 77(1*) + \ J R \ogH sui 2)(t, e ix , e")da * da{x). 

Here, Hsu (2) is the central heat kernel of SU(2) and e %x is short for the diagonal ma- 
trix D{e LX 1 e~ lx ). The element a is the principal element of type p of SU(2), namely 
a = D(e 2 , e 2 J. 

As discussed in jF] (Proposition 1.3), the central heat kernel of 577(2) at these special 
values is given by 

where 



^(2)(t,e fa ,eT) = o ^ M J^ (6) 



0i(z, t) = 2 ^(-l) n sin{(2n + l)z}e-^ n+l/2)H (7) 

n=0 

is Jacobi's theta function. MacDonald's identities are more often stated in terms of Jacobi 
theta functions, but we find it easier to work directly with heat kernels. We note that the 
factor of 7r in the exponent 7r(n + l/2) 2 t is responsible for the appearance of nt in many 
expressions to follow in the heat kernel. 

At first sight, this result seems to explain the normalization -™ logi^y. However, the large 
N limit conjectures concern the scaling limit of the central heat kernel under t — > This 
puts into play a simultaneous limit process in da^*daN —>■ da* da and in the asymptotics of 
theta functions. We find that for our cases of the problem, this rescaling changes the growth 
rate of the free energy. 

Theorem 1.3. Suppose that e kN is a sequence such that da^ — ► da. Then, 
F N (a N ,e k »\^) = \ J sl \ogH sui2) (A/2N,e ix ,ef)da N *da^(e ix ) - >gr/(§) + 0(1/ N) 

~ -f {f gl J min{d{e ix , e in / 2 ), d{e ix , e~™/ 2 )} 2 da * da(e ix ) - 

where d(e lx ,e w ) is the distance along S 1 . 

The minimal distance above arises as the distance d(C(a),e tx ) in SU(2) of the diagonal 
element D(e lx , e~ rx ) to the conjugacy class C(a) of a of the principal element of type p, which 
is conjugate to D(e t7T ^ 2 , e~ m l 2 ). 
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We note that there are two terms of opposite sign in the leading order term. We note that 
the terms cancel when U-i = a = Uo, since then da * da = ^ and the first term reduces to 

1 z ' 2-7T 

— yH^:^!(i) 3 ~~ T^J" = ®' This a PP ears to be the kind of case studied in |GMlj . But the 
two terms cannot cancel in all cases, and indeed, the leading term does not not cancel in the 
simplest case, where = for all N, i.e. where U 2 = Id. We then have (see |Flj or |Fj . 
Proposition 1.2) 

H su{N) (t, a N , I) = e -*»suW tfiu^snsum (8) 

where rj(t) is Dedekind's 77-function. As is easy to see (and will be verified below), the 
eigenvalue distribution of tends to while that of / is obviously 5\. In this case, the 
asymptotic mass equals 1 and the continuous term equals zero. We separate out this special 
case since the result is most easily checked on this example: 

Proposition 1.4. When U\ = a N and U 2 = I, so that o\ = d0, a 2 = 5\, then 

A- 2 \ogZ su{N) (e^, 1\A) = -f {£} - | log(w) - A / 48N + °( e ~ cN )- 
Note that Theorem 11.31 reduces to Proposition II .41 when da = 5q. 

Some final comments on the anomalous extra factor of N. The prediction that log Zn (A) 
should have a limit determined by a variational problem is one of many predictions of this 
kind in field theory and statistical mechanics and it seems very strange that an extra factor 
of N should appear in our calculations. It is a mystery how it would appear in the graphical 
or diagrammatic calculations which initially suggested that the large N limit of gauge theory 
is a string theory, i.e. without using (JTJ). It seems best to leave it to the string theorists to 
decide how much the anomalous factor of iV affects the picture of the large N limit of gauge 
theory as a string theory. 

We are on somewhat safer ground in trying to account for the extra factor of iV in the 
formal calculations based on (fl|). From discussions with M. Douglas and V. Kazakov, it seems 
plausible that the anomaly is due to the unusual behavior of the character values Xni ) a t 
the principal element of type p. It was proved by Kostant |Koj that the only character values 
at this element are Xn{ a ) = ~ 1, 0, 1 and it is reasonable to expect that the value oscillates 
regularly between these values. It appears that the physics predictions implicitly assumed a 
less oscillatory behaviour in character values, and in particular less oscillation in the signs of 
character values. The heuristic calculations in |GM1[ IGM2| IKWj of the large N limit of the 
free energy on the cylinder were based on special cases (such as U 2 = U*) which do not have 
such oscillations in sign. The sign oscillation causes much more cancellation than expected, 
and this could explain why our asymptotics have the form e _Ar3 " rather than e _Ar2= . These 
special values Xi?( a ) might also be inconsistent with Matytsin's character asymptotics, and 
we plan to check this in the future. 

This paper is organized as follows. In Section we review MacDonald's identity and 
associated objects which we will use in the proof of the main result. The main results are 
proved in Section El In Section we discuss some aspects of the proof and mention some 
other results which could be proved by the methods of this paper. 

Acknowledgements The author would like to thank M. Douglas, P. Etingof, V. Kazakov, 
B. Kostant, N. Reshtikhin, T. Tate, and C. Vafa for going over some of the details of the 
calculations and for comments (quoted above) on the relations between our results and the 
predictions in the physics papers. 
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2. Background on MacDonald's identities 

In this section, we review the rj function and MacDonald's identity. Our main references 
are the articles [F?] of I- Frenkel and those jFJ IFlj of H. Fegan. Since the articles employ 
very different, possibly confusing, notational conventions, we highlight the main ones: 

• In |Frj . the 77-function is defined in a non-standard way in the lower half-plane (see 
§4.4) and consequently the heat kernel is evaluated at time t = Airib where 96 < 0. 
We will quote results of jFrj in terms of t and we will use the usual definition of r\ as 
a modular form on the upper half plane. 

• In jFfj, the special element of type p is denoted e 4mp , but the exponent refers to the 
dual element of the Cartan subalgebra rather than the linear functional p. To avoid 
confusion, we denote it simply by a as in jFj. 

• Fegan works with the Schrodinger equation rather |Flj (1.1) rather than the heat 
equation, so our formulae differ from his in that his it is our t. 



2.1. i] function. The Dedekind eta-function jUE] is defined by 

v {z) = e 2mz/24 n£° =1 (l - e 2 ™ 2 ), Qz > 0. (9) 
It is a modular cusp form of weight 1/2, i.e. it satisfies 

V ( 7 z) = e^izY^z), if 7 e SL{2, Z), 

where #(7) is a certain multiplier which we will not need to know in detail (see [Q, §2.8). 
When jz = — 1/z, we have 

v(-~) = (iz) 1/2 v(*)- 

z 

In studying the free energy, we are particularly interested in logrj(iy) as y — * 0+. The 
transformation law for logr](iy) goes as follows |7Q: For real numbers y > we have: 

\ogr](iy) = - ^ - \ log y + YZ=i £ 1=^75 • (10) 

2.1.1. Jacobi's theta function. Although we mainly use the heat kernel parametrix on SU (2) 
to obtain asymptotic formulae, the alternative in terms of Jacobi's theta function (jUJ) can be 
used to check the details. The small time expansion of the heat kernel can be obtained form 
Jacobi's imaginary transformation law: 

9i(-x/t,-1/t) = i(-iT) 1/2 e- lx2/T 6(x,r). 



Thus, we have: 



Here, we use that 



iA 2N _^a_ _ n 2N 

- — ----- - jr - 1— =>• e 2iv _, e 



The large N or small time asymptotics now follow from (JJJ). 
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2.2. Central heat kernels. Let G be any compact, connected Lie group. We denote by R 
the root system of (gc, he), where h is its Cartan subalgebra. We denote by R + the positive 
roots, by P the lattice of weights and by P ++ C P the dominant weights. 

We recall that the eigenvalue of the Casimir operator (bi-invariant Laplacian) A of G in 
the representation with highest weight A is equal to 

A|v A = (||A + p|| 2 -||p|| 2 )Jrf| yA , p=l/2 «■ 

The fundamental solution of the heat equation is given by k(t, x, y) = v(xy~ x , t) where 

v(g,t) = £ (dim^XA^e-^ll^l'-ll""'), (12) 
AeP++ 

As in jEj, §4.3 it then follows that the central heat kernel is given by 

H G (t,h,k)= XA(e' l )xA(e- fe )e-*/ 2 (ll^ll 2 -ll^ 2 ). (13) 
\eP++ 

This is simply a different notation for (J2J). 

2.3. SU(2). As was recognized clearly by H. Fegan and others, the central heat kernel of 
SU(2) plays an important role in MacDonald's identities. There is one positive root a, 
which can be identified with 1 if we choose it as the basis of the Cartan dual subalgebra. 
Then p — \ and the Killing form is B(x,y) = \xy. In this case, the principal element a 

of type p has eigenvalues e ±1 ^ . The weight lattice is |Z. The character of the irreducible 

representation of highest weight A is X\( x ) — _ We have: 

— 1, A is an odd integer 
X\(a) = { 0, is not an integer 
1, is an even integer 

The eigenvalues of the Casimir are c(A) = |A(A + 1). 

The central heat kernel H (t, a, y) at the special point a may be expressed in terms of 
Jacobi's theta function as in ©• As a special case of MacDonald s identities, we further 
have (cf. [EE], (3.12)): 

H su{2) (t,a,l) = (e-^ 12 V (it)) 3 . (14) 

2.4. MacDonald identities. We briefly review MacDonald's identities for a compact, semi- 
simple, simply connected Lie group G and its relation (proved by Kostant |Koj ) to the central 
heat kernel. They are thus valid for SU(N). We follow |F1 IFr]. 

MacDonald's identity involves a conjugacy class C a of special elements of G, which we 
will denote by a, which are conjugate to e 2Xp where x p is dual under the Killing form to 
p. Such elements are called 'principal elements of type p'. A principal element of type p 
is characterized in iKoj (see 1.3 p. 181) as a regular element such that the order of Ad(a) 
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equals h (the order of the Coxeter element). One has h = N for SU(N), so the eigenvalues 
of Ad(a) must be distinct iVth roots of unity. Hence 

OjV -> r-, (15) 

i.e. the eigenvalues of become uniformly distributed in the large N limit. A detailed 
description of such elements can be found in |Ko[ IF1[ lAFj . 

In Kostant's formulation, as described by Frenkel, the MacDonald's identities take the 
form ([El, Proposition (4.4.5)) 

e(k,4mt) = *(k) XA(a)XA(e- fc )e- t/2||A+ ^ 112 , (16) 
AeP++ 

where (see |Prj . (1.1.8) 

a(k)=U aeR+ (e^ 2 -e-^ 2 ), 
where 9(k,it) is the theta-function defined in Definition (4.4.1) of [Frj: For t > 0, 



;i7) 



9(k,it) = e- 2 -IHI 2 n aei?+ (e^-e-^) 

• n~ =1 (i - e- 2nnt yu aeR (i - e - 2 ™'+< a > fc >) 

For our purposes the key formula is the following ((El, Proposition (4.4.4)): 

9(k,it) =r ] (it)-\ R +\ +e Tl a€R+ 6 1 ((a,k)),tt). (18) 

We put 9(x, it) = 2e-t/i'sinTTx ' ^ e following is the version of MacDonald's identities proved 
in Theorem 1.5 of |Flj (together with its Proposition 1.3). As above, an element e ld is 
identified with a diagonal element D{e l9 1 e~ td ) of SU(2). 

Lemma 2.1. Let G be compact, connected, semi- simple and simply connected, and let a be 
an element of type p. Then 

H G (t,a,e~ k ) = {e^ V m- mi "*^Hsup)(t,^\e*^ 
Proof. Puttting together (JTfiJl . (fTHj) and (|T3^1. we get 

H G (t,a,e- k ) = a(k)~ 1 e A/4N ^ 2 9(k, it) 

(19) 

= e t IWI a ^(i*)- |Afl+£ n a6At 0((a,ife),zt). 
The stated result now follows from (0). □ 

As mentioned above, the simplest case (j2"Uj) comes from combining Lemma f2. II and ([14)): 

H G {t,a,I) = e- dimGt / 24 r](it) dimG , (20) 
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3. Large iV limit of Z S u(N){a,U\A): Proof of Theorem 11.31 

3.1. Simplest case. We first consider the simplest case ()2Uj) . Since k = (i.e. XJ<i = I), we 
are essentially dealing with YM 2 on the disc. This case is discussed in detail in |GM2j . §4. 

Proposition 3.1. When U\ = a and U 2 — l, so that o\ = d6, a 2 =81, then 

^\ogZ su{N) (a,l;A) = -f £ - \ log(^) - A/A8N + 0(e~ cN ). 

Thus, no phase transition occurs. 

Proof. Macdonald's identity gives: 

Z su(N) (a,l;A) = e- dim5f/ W A / 247V r ? (zA/2Ar) dim5C/ W, . (21) 

The free energy is then 

±- 2 \ogH su(N) (A/2N,a,I) = dimS ^ { _^_ + teg v {iA/2N)}. (22) 

We note that dimS N u 2 {N) = \ + 0(±). We substitute y = A/2N in the right side of (HH to 
get: 

/ . / n 2N TT 1 , , A . ^ 1 1 

log,(M/2A0 = - -log(-) + £ - t.^wx < 23 > 

m=l 

The summand of the sum J2m=i m i- *n™/a 1S smaller than \/2e~ AN ^ K ' A ^ m ) for N suffi- 
ciently large, so we may bound the sum by the geometric series and obtain a bound of 

- — } m */A) ~ 1 < Ce~ mT/A \ for N sufficiently large. 
Thus, the sum is an exponentially small correction, and we have 

^\ogZ su(N) (a,l-A) = ^^{-A/48iV-|f -Ilog(^)} + 0(e-^). 



(24) 

\ — A /ARM 

24 A 2 1W &V2A^ 



T2V -|M4)-^/48iV. 



□ 



We see clearly the anamolous factor of N. Also, we see that there is no phase transition, 
despite the fact that the two eigenvalues densities, 5(1), dd have very different support prop- 
erties. In almost exactly this case, the lack of phase transition is predicted by Gross-Matytsin 
in |GM2j . §4. They explain that when the support of <j\ is the whole circle and the support 
of o~ 2 is a single point, then the system is always in the strong coupling phase. 

They further explain that when o~\ = 5(6 = 0), the solution of the Hopf equation can be 
obtained as the solution of the integral equation 



A' J v ' ' A J 6-6' -tf(t,6) 
In our case, 01 = 1 so the equation is simply 

(l-j)f(t,6) = ^--\og(6-tf(t,6)). 
It would be interesting to solve this equation and check that the solution is analytic in A. 
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3.2. Proof of Theorem (11. 2 j) . We may take e k to be a diagonal matrix with entries 
e 27,i\j(N)_ p irgt) we put t _ Aj2N. By Lemma Owe have: 

^\ogH su{N) (t,a,e k ) = ^{(-\R + \ + N)[-^+log V (-it/A7r)} 

(25) 

+ wE a> o l ogH su{2) (t,e^,e^) 

We note that both sides of this equation are real, so that we take take the real part 3? without 
changing the equation. 

We now specialize to the case of SU (N) . Its roots are ej — ej and its positive roots satisfy 
i < j. Hence, {a, k) = A« — Xj. Hence 

£ aei?+ hgH su{2) (t, e*<«.*>, e-/2) = ^ hgH su{2) {t, e^~H e^ 2 ) 

(26) 

= iE^log^c/(2)(t,e^- A ^,e-/ 2 ). 

In the last equality we use that H SU (2)(t, e i{x ^^\ e i7r/2 ) = H SU (2)(t, e i(A ^ Al) , e i7r/2 ). 
We further have 

logH su{2) (t, e^"H e™' 2 ) = J2 logH su{2) (t, e<^\ e^ 2 ) - N \ogH su{2) (t, 1, e^ 2 ). 

Since £at ~ iV for SU(N), we have 
W T, a&R+ logH su(2) (t, e^ k \ e-/ 2 ) = J sl J sl \ogH SU (2)(t, e ix , e^ 2 )da N (e^-^)d^(e lx ') 



^log# sc/(2) (t,l,e-/ 2 ) 



where 



If 



then 



J sl log H su(2) (t, e ix , e™' 2 )da N * da N (e" 
-±logH su{2) (t,l,e^ 2 ) 

da N *da^(e ix ) = [ da N {e i{x - x>) ) da N (e ix '). 
's 1 



In 



da^ * (icrTv ^ da * da, 

since the Fourier coefficients of the left side tend to those of the right side. Thus, we obtain 
the stated limit. 

3.3. Proof of Theorem 11.31 We now re-do the calculation but make the scaling t = 
We thus have 

F su{N) {a,e k ;A) = \ J gl log H SU (2)(^, e ix , e i7r/2 )da N * da^{e ix ) 

(27) 

-±logH su{2) (^,l,e^ 2 ) 
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We now obtain the limit by using the uniform off-diagonal asymptotics of the central heat 
kernel. 

The central heat kernel is related to the actual heat kernel by 

log H SU (2)(t,e tx ,a) = log / k SU (2)(t,e lx ,g~ 1 ag)dg 

JSU{2) 

where ksu(2){t,x,y) is the heat kernel. Therefore, we are interested in the uniform asymp- 
totics of 

A f A 

log H su{2) {—, e ix , a) = log / k su{2 ){^, e ix , g" 1 ag)dg (28) 

J SU (2) 

in x for each A. There exists a uniform heat kernel parametrix for ksu{2) given by: 

ksu(2)(t, u, v) ~ t- 3 / 2 e~^^V(t, u, v) (29) 

where V(t, u, v) ~ YlJLo v)P. The amplitude V is to leading order the volume density 
in normal coordinates. We will never be evaluating it at a pair of conjugate points, so it is 
uniformly bounded below by a positive constant and introduces only lower order terms into 
the logarithm. See for instance jKaj for general results of this kind. Instead of a parametrix, 
the reader might prefer to use the (JBJ) in terms of Jacobi's theta function. This formula and 
also (jl4j) can be used to check the various constants in the formula. 
Thus, we have 

log H su(2) (^,^,a) ~log{(^) 3 / 2 I e-^ d(elx ^ 2 V(A/2N,e^,g- 1 ag)dg}, (30) 

where V(A/2N, e ix , g^ag) is a semiclassical amplitude in N. The asymptotics are deter- 
mined by the minimum point of the phase d(e lx , g _1 ag) 2 , namely by the distance d(e lx ,C a ) 
from e lx to the conjugacy class of a. We note that the conjugacy class C(a) = {g _1 ag : g G 
SU{2)} is a great (equatorial) 2-sphere of radius 7r/2 from (the north pole) /. 

There is a somewhat different expansion accordingly as e lx = ±1, e lx G C(a) or for e lx not 
of this form, which we will call a general element. If e lx is a general element, then there is a 
unique closest point g~ x ag and the phase is non-degenerate, so we obtain 

Isu { 2) e- 2 ^ x ^ la ^V(A/2N, x, g- x ag)dg 

(31) 

~ A' 3 / 2 e -f^(e«C a )3 . 1 

^/det Hess{A~ 1 d(e lx ,g~ 1 ag) 2 ) 

We note that the powers of A cancel. If e lx G C(a), then since e lx also represents a point in 
the maximal torus, e tx = a or e lx = a -1 , hence the unique point of minimal distance to e lx 
in C (a) is of course e lx itself. There is no essential change in the calculation except that the 
exponent vanishes. We now consider the behavior of -, 1 . = as e ix — > ±1. In 

W det Hess(d(e lx ,g~ 1 ag)' 2 ) 

the case e lx 1), for instance, det Hess{A~ 1 d(e lx , g~ 1 ag) 2 ) ~ |x| 3 and hence 

log det Hess(d(e lx , g~ l ag) 2 ) ~ log \x\ as x — > 0. 

At e lx = ±1, the entire C(a) becomes a critical manifold for the phase. Exactly at the poles, 
we have 

H sm (t,±l,a) = k sm (^,±l,a) ~ (^3/2^2^/4^. (32) 
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Thus, as e ix — ► 1, 



.4 



\ x \-V2 e -^d(e-,C a )^ x o 



^(^.e^a)- (33) 

- A 1 /Af^3/2 e -27V7r 2 /47rA x = 

Similarly at x = tt. We note that the exponent is continuous and only the power of N changes 
at the special points x = 0, tt. 

Since we are interested in log asymptotics, the factor e~^ d ^ x ' Ca ^ is dominant as long 
as the remainder can be integrated against da^ * da^. If this measure has a point mass at 
either or tt, there are singularities at x = 0, tt in the coefficients of the asymptotics which 
are not integrable da * da. We now discuss this point in detail. 

We fix a positive constant C > and break up the integral (|33j) as 



Jsi\i-c/n,c/n\ui*-c/n,*+c/N] l °g H su(2){^, e %x , e m l 2 )da N * da N (e lQ 

+ I[-C/N,C/NMn-C/N,7r+C/N] lo S H su{2) ( ^, e™ , e m ^)da N * da N (e lc 
It follows from the pointwise asymptotics in (jHSJ) that 

Is i \[-c/N,c/N]u[ ?r -c/N^+c/N\ l °g H su(2)(4fi e ™ > e**/ 2 )da N * da N (e l 

~ J^\[-c/iv,c/iV]u[ 7 r-c/iV l7 r+c/iv] e") 2 ^ * ^(e^). 



(34) 



(35) 



'iV' AN 

asymptotics of the heat kernel rather than its pointwise asymptotics. Since the details are 



For the second integral of over [— £ ^] and [7r — C /N,tt + C /N], we use the scaling 

an its 

similar for both intervals we only carry them out around the interval [— §;]. Namely, we 



write x = jfe with \u\ < C. We have: 



d{e^ l9 - l agf = (^) 2 + ^Q(u, N, A, g^ag) 



for a bounded smooth function Q, hence 

ISU(2) 



log H su{2) (^,e^,a) ~ {(f) 3 / 2 e-^/ 5a(2) e-^<^^ 



~ (f)3/2 e -^4[ g ( M ,A) + ig 1 ( M ,A) + ---) 

(36) 

where q = e - Q ( u ' A '°) is strictly positive. Substituting (jHHjl into the second term of (J33|) . we 

get 



7rA 4 



3 N f u f u 

2 lo g(^[)] J c da N * da N (—) - J Q(u,A)d<r N * da N ( — ). (37) 

Since Q is bounded and daN(j^) * da^ijf) has mass at most 1, the last term is 0(1) as 
N — > oo. Thus, the first term of (}3*Tj) dominates. Thus, we have 

Il-c/N,c/N]uin-c/N,n+c/N] lo S H sum(m> e " "> e^ 2 )^ * da^(e ix ) 

(38) 



MACDONALD'S IDENTITIES AND THE LARGE N LIMIT 



OF YM 2 ON THE CYLINDER 13 



We note that d(C a , 1) = d(C a , — 1) = |. Since d(C a ,e lx ) is continuous, we may rewrite 
fj37|) - (j38j) to leading order as 

2N r c ' N 



7TA 

We add this back to (J35J) to obtain the leading order term 

2N 



d(e ix ,C a ) 2 da N *da N . (39) 

-C/N 



TlA 



d(C(a), e lx yda N * da N (e lx ) (40) 



plus the canonical terms log H S u(2) (A/2 N, 1, a) and logr/(|^) which are independent of 
do - . We then recognize that d(e lx ,C a ) = mm{d(e lx , e ±ln / 2 )}, completing the proof. 

□ 

Remark: 

• It would be interesting to know when the leading order term cancels, but even this 
cancellation would not cure the anomoly, since there is no control over the rate of 
the weak convergence da^ — > da. Hence, there is no well-defined growth rate of the 
lower order terms. 

• As noted above, the exponent in ()33|) is continuous. Since we are taking the logarithm, 
only the exponent is important to leading order, and that is why we obtain a unified 
formula, even when da has a point mass at e lx = ±1. 



4. Final Remarks 

Some final comments and remarks. 

• It would of course be desirable to remove the assumption that one of the matrices Ui 
should be a prinicpal element of type p. In |F2j . H. Fegan states a product formula 
for the full fundamental solution, but unfortunately this formula is erroneous. It 
would be interesting to see if there exists some form of MacDonald's identities for 
the full central heat kernel. 

• It would also be desirable to extend the results of this paper to more general groups 
such as U(N), SO(N), Sp(N). We specialize to = SU(N) because it amply illus- 
trates the main point of this paper. Further, MacDonald's identities pertain to simply 
connected compact semi-simple groups, and presumably require some modifications 
for U(N),SO(N) etc. 

• If we ignore the unexpected growth in N, the phase transition conjecture suggests 
that the limit free energy should equal 

S(d9, da\A) — | J sl J sl log | sin ^■\da(6)da((f)) — | J gl J sl log | sin ^-\d6d(j), 

The last two terms can be obtained from our limit formula by unravelling the denom- 
inator in the logarithm. Thus, aside from orders of magnitude, the phase transition 
conjecture appears to state in these cases that 

S(d9, da\A) = ^-r [ d{C{a), e ix fda * da~(e ix ). 
2A J s i 
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We have not been able to compare our results directly with the physics predictions 
by solving the complex Burghers equation. 

• In addition to considering 'pointwise' asymptotics of the partition function, one could 
consider weak asymptotics where one averages in various ways over the variables 
U\,U2, considers the variance of the integrals and so on. 

• As mentioned in the introduction, it is possible that an anomaly also occurs in 
Matytsin's conjectured asymptotics |M] at the special character values Xr(°)- 
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